home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / modelers / geomview / source.lha / Geomview / src / lib / gprim / geom / iterate.c < prev    next >
Encoding:
C/C++ Source or Header  |  1993-04-18  |  6.2 KB  |  288 lines
  1. /* Copyright (c) 1992 The Geometry Center; University of Minnesota
  2.    1300 South Second Street;  Minneapolis, MN  55454, USA;
  3.    
  4. This file is part of geomview/OOGL. geomview/OOGL is free software;
  5. you can redistribute it and/or modify it only under the terms given in
  6. the file COPYING, which you should have received along with this file.
  7. This and other related software may be obtained via anonymous ftp from
  8. geom.umn.edu; email: software@geom.umn.edu. */
  9.  
  10. /* Authors: Charlie Gunn, Stuart Levy, Tamara Munzner, Mark Phillips */
  11.  
  12. /*
  13.  * iterate.c -- iteration over OOGL hierarchies
  14.  */
  15.  
  16. #include <geomclass.h>
  17. #include "listP.h"
  18. #include "tlistP.h"
  19. #include "instP.h"
  20. #include "sphereP.h"
  21. #include "discgrpP.h"
  22.  
  23. #define    UNKNOWN    0
  24. #define    LIST    1
  25. #define    TLIST    2
  26. #define    INST    3
  27. #define SPHERE    4
  28. #define DISCGRP 5
  29. #define OTHER    6
  30.  
  31. #define    VALID    0x13ac2480    /* Magic number */
  32.  
  33. struct istate {            /* Iteration state structure: per nest lvl */
  34.     struct istate *parent;
  35.     Geom    *g;
  36.     int    kind;
  37.     int    seq;        /* Seq number within Tlist */
  38.     Transform Ti;
  39. };
  40.  
  41. struct GeomIter {
  42.     struct istate *stack;
  43.     int flags;
  44. };
  45.  
  46. static struct GeomIter *itpool;
  47. static struct istate *ispool;
  48.  
  49. #define    NewIter(it)    if(itpool) { it = itpool; itpool = NULL;} \
  50.             else it = (GeomIter *)malloc(sizeof(GeomIter))
  51. #define    NewIstate(is)    if(ispool) { is = ispool; ispool = ispool->parent; } \
  52.             else is = (struct istate *)malloc(sizeof(struct istate))
  53.  
  54. #define    FreeIter(it)    if(itpool) GeomFree(it); else itpool = it
  55. #define    FreeIstate(is)    (is)->parent = ispool, ispool = is
  56.  
  57. GeomIter *
  58. GeomIterate(g, flags)
  59.     Geom *g;
  60.     int flags;
  61. {
  62.     GeomIter *it;
  63.     register struct istate *is;
  64.  
  65.     NewIter(it);
  66.     it->flags = (flags & 15) | VALID;
  67.     NewIstate(is);
  68.     it->stack = is;
  69.     is->kind = UNKNOWN;
  70.     is->g = g;
  71.     is->parent = NULL;
  72.     return it;
  73. }
  74.  
  75. /*
  76.  * Iterate over an OOGL hierarchy, producing one Transform per call.
  77.  * Returns 1 normally, 0 when no more Transforms remain.
  78.  * Destroys the GeomIter upon returning 0.
  79.  *
  80.  * To make this work as a coroutine we effectively unroll the recursive
  81.  * hierarchy traversal, so please excuse the goto's.
  82.  *
  83.  * For speed, we act as a friend to the List, TList, Inst, and Sphere
  84.  * classes.   XXX
  85.  */
  86. NextTransform(it, T)
  87.     GeomIter *it;
  88.     Transform T;
  89. {
  90.     register struct istate *is;
  91.     Geom *g;
  92.  
  93.     if(it == NULL)
  94.     return 0;
  95.     if((is = it->stack) == NULL)
  96.     goto gone;
  97.  
  98.     g = is->g;
  99.  
  100.   again:
  101.     switch(is->kind) {
  102.     case UNKNOWN:                /* Determine type */
  103.   discover:
  104.     if(g->Class == TlistClass) {
  105.         if(g->magic != TLISTMAGIC) {
  106.         GeomError(0,"NextTransform: Surprise Tlist %x", g);
  107.         goto pop;
  108.         }
  109.         is->kind = TLIST;
  110.         is->seq = 0;
  111.         goto isTLIST;
  112.     }
  113.     if(g->Class == DiscGrpClass) {
  114.         if(g->magic != DISCGRPMAGIC) {
  115.         GeomError(0,"NextTransform: Surprise DiscGrp %x", g);
  116.         goto pop;
  117.         }
  118.         is->kind = DISCGRP;
  119.         is->seq = 0;
  120.         goto isDISCGRP;
  121.     }
  122.     if(g->Class == ListClass) {
  123.         if(g->magic != LISTMAGIC) {
  124.         GeomError(0,"NextTransform: Surprise List %x", g);
  125.         goto pop;
  126.         }
  127.         is->kind = LIST;
  128.         goto isLIST;
  129.     }
  130.     if(g->Class == InstClass || g->Class == SphereClass) {
  131.         register Inst *inst = (Inst *)g;
  132.  
  133.         if(inst->magic != INSTMAGIC && inst->magic != SPHEREMAGIC) {
  134.         GeomError(0,"NextTransform: Surprise Inst %x", g);
  135.         goto pop;
  136.         }
  137.         if(inst->tlist == NULL) {
  138.         /*
  139.          * Just use this transform.  We're a leaf, so return it now.
  140.          * Also pop stack for next iteration.
  141.          */
  142.         if(is = is->parent)
  143.             TmConcat(inst->axis, is->Ti, T);
  144.         else
  145.             TmCopy(inst->axis, T);
  146.         FreeIstate(it->stack);
  147.         it->stack = is;
  148.         return 1;
  149.         } else {
  150.         /*
  151.          * Inst has a child.
  152.          * The axis field is *ignored* in this case; just
  153.          * traverse the child.
  154.          */
  155.         is->g = g = inst->tlist;
  156.         /*is->kind = UNKNOWN;*/
  157.         goto discover;
  158.         }
  159.     }
  160.     /* Other objects aren't relevant for NextTransform(), so discard. */
  161.     goto pop;
  162.  
  163.     case TLIST:
  164.       isTLIST:
  165.       {
  166.     register Tlist *TL = (Tlist *)g;
  167.     register float (*Tp)[4];    /* ought to be ``Transform *Tp'' */
  168.  
  169.     if(is->seq >= TL->nelements)
  170.         goto pop;
  171.  
  172.     Tp = TL->elements[is->seq++];
  173.  
  174.     if(TL->tlist != NULL && !(it->flags & SHALLOW)) {
  175.         if(is->parent)
  176.         TmConcat(*Tp, is->parent->Ti, is->Ti);
  177.         else
  178.         TmCopy(*Tp, is->Ti);
  179.         g = TL->tlist;
  180.         goto push;
  181.     }
  182.  
  183.     /* We're a leaf -- return a matrix now */
  184.     if(is->parent)
  185.         TmConcat(*Tp, is->parent->Ti, T);
  186.     else
  187.         TmCopy(*Tp, T);
  188.     return 1;
  189.       }
  190.  
  191.     case DISCGRP:
  192.       isDISCGRP:
  193.       {
  194.     register DiscGrp *dg = (DiscGrp *)g;
  195.     register float (*Tp)[4];    /* ought to be ``Transform *Tp'' */
  196.  
  197.     if(is->seq >= dg->big_list->num_el)
  198.         goto pop;
  199.  
  200.     Tp = dg->big_list->el_list[is->seq++].tform;
  201.  
  202.     /* if it's another discrete group ... */
  203.     if(dg->geom != NULL &&
  204.         dg->geom->magic == DISCGRPMAGIC && !(it->flags & SHALLOW)) {
  205.         if(is->parent)
  206.         TmConcat(*Tp, is->parent->Ti, is->Ti);
  207.         else
  208.         TmCopy(*Tp, is->Ti);
  209.         g = dg->geom;
  210.         goto push;
  211.     }
  212.  
  213.     /* We're a leaf -- return a matrix now */
  214.     if(is->parent)
  215.         TmConcat(*Tp, is->parent->Ti, T);
  216.     else
  217.         TmCopy(*Tp, T);
  218.     return 1;
  219.       }
  220.  
  221.     case LIST:
  222.      isLIST:
  223.       {
  224.     /* Might need some intra-List state too if List format
  225.      * gets fancier -- if we have one List object which contains a
  226.      * bundle of Geom's.
  227.      */
  228.     register List *L = (List *)g;
  229.  
  230.     if(L == NULL)
  231.         goto pop;
  232.     while(L->car == NULL) {
  233.         L = L->cdr;
  234.         if(L == NULL)
  235.         goto pop;
  236.     }
  237.     g = L->car;
  238.     is->g = (Geom *)L->cdr;
  239.     if(is->parent)
  240.         TmCopy(is->parent->Ti, is->Ti);
  241.     else
  242.         TmIdentity(is->Ti);
  243.     goto push;
  244.       }
  245.  
  246.     default:
  247.     GeomError(1,"NextTransform: called with invalid GeomIter (%x)", it);
  248.     return 0;
  249.     }
  250.  
  251.   push:            /* Push g as a new element */
  252.             /* is->Ti *must* be set correctly before getting here! */
  253.     NewIstate(is);
  254.     is->g = g;
  255.     is->parent = it->stack;
  256.     it->stack = is;
  257.     goto discover;    /* Recurse */
  258.  
  259.   pop:
  260.     is = is->parent;
  261.     FreeIstate(it->stack);
  262.     it->stack = is;
  263.     if(is != NULL) {
  264.     g = is->g;
  265.     goto again;    /* Recurse -- type already determined */
  266.     }
  267.     /* stack empty -- fall into gone: */
  268.  
  269.   gone:
  270.     it->flags = 0;
  271.     it->stack = NULL;
  272.     FreeIter(it);
  273.     return 0;
  274. }
  275.  
  276. void
  277. DestroyIter(it)
  278.     GeomIter *it;
  279. {
  280.     if((it->flags & 0xfffffff0) == VALID) {
  281.     it->flags = 0;
  282.     it->stack = NULL;
  283.     FreeIter(it);
  284.     } else {
  285.     GeomError(1,"DestroyIter -- already destroyed %x", it);
  286.     }
  287. }
  288.